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Abstract 

An effective monopole action is derived from vacuum configurations after 

abelian projection in the maximally abelian gauge in SU(2) QCD. Entropy 

dominance over energy of monopole loops is seen on the renormalized lattice 

with the spacing b > b c ~ 5.2 x 10 _3 A^ 1 when the physical volume of the 

system is large enough. QCD confinement may be interpreted as the (dual) 

Meissner effect due to the monopole condensation. 
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To understand confinement mechanism is very important but still unresolved problem in 
particle physics. A promising idea is that the (dual) Meissner effect due to condensation of 
some magnetic quantity is the color confinement mechanism in QCD 0,0]. This picture is 
realized in the confinement phase of lattice compact QED [Q-|5| . Especially interesting are 
the following facts. 

1) A dual transformation can be done, leading us to an action describing a monopole 
Coulomb gas. Monopole condensation is shown to occur in the confinement phase 
from energy-entropy balance. 

2) The monopole contribution alone can reproduce the full value of the string tension 

In the case of QCD, there is a difficult problem. We have to find a color magnetic 
quantity in QCD. In this respect, the 'tHooft idea of abelian projection of QCD |L0| is 
very interesting. The abelian projection of QCD is to perform a partial gauge-fixing such 
that the maximal abelian torus group remains unbroken. Then QCD can be regarded as a 
U(l) x U(l) abelian gauge theory with magnetic monopoles and electric charges, 't Hooft 
conjectured that the condensation of the abelian monopoles is the confinement mechanism 
in QCD []10[ . There are, however, infinite ways of extracting such an abelian theory out of 
SU(3) QCD. It seems important to find a good gauge in which the conjecture is seen to be 
realized. 

A gauge called maximally abelian (MA) gauge has been shown to be very interesting 



11-14 . In the MA gauge, there are phenomena which may be called abelian dominance 



j~2|,?l. Moreover the monopole current k^s) can be defined similarly as in compact QED 
. It is shown in the MA gauge that the abelian monopoles are dense and dynamical in 



the confinement phase, whereas they are dilute and static in the deconfinement phase |TJ| . 

The abelian dominance suggests that a set of U(l) invariant operators are enough to 
describe confinement in the MA gauge. Then there must exist an effective U(l) action 
describing confinement. Let us try to derive an effective U(l) action on the dual lattice as 
done in compact QED from monopole distributions given by Monte-Carlo simulations. 



Swendsen jHJ developed a method of determining an action from a given ensemble of 



configurations. The method can be applied to our case. A theory of monopole loops is given 
in general by the following partition function 

oo 

2 = 01 E )(n<W),o)exp(-£^M), 

*>M fc M (s)=-oo s i 

where <9' is a backward derivative on a dual lattice and fc^(s) is the conserved integer- valued 
monopole current, /j is a coupling constant of an interaction Si[k]. Since the dynamical 
variables here are k^s) satisfying the conservation rule, it is necessary to extend the original 
Swendsen method by considering a plaquette (s', fjf, v') instead of a link. Introducing a new 
set of coupling constants {fi}, define 

g _ Em=-co £fjjy ex P(~ Si /ffiffij) ^ 
exp(-Ei/ J 5'i[ fc/ ]) 
where fe^(s) = /c M (s) + M{8 a>s i8 lllJ ; + £ SjS , 

-i-p,' "liv' "SjS'+u' $mj,' - 5 S , S '^')- When all /j are equal 
to fi, one can prove an equality (Si) = (Si), where the expectation values are taken over 
the above original action with the coupling constants {fi}- Otherwise, one may expand the 
difference as follows: 

(St ~ Si) = ^(SWj - SS ,){]) - fj), (2) 
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where only the first order terms are written down. This allows an iteration scheme for 
determination of the unknown constants fi. 



Practically we have to restrict the number of interaction terms [I6|. We adopted 12 types 
of quadratic interactions listed in Table | in most of these studies [T7J . 

First we applied the method to the Villain form of U(l) lattice gauge theory [ 118 , 19 1 , since 
the Villain model is reformulated exactly as a theory of monopole loops having an action 

J2fiSi[k] = 2tt 2 A, £ k,(s)D(s - s')k,(s'), (3) 

i s,s',fj, 

where D(s) is the lattice Coulomb Green function [[|. 

Using the Villain action, we generated 100 gauge field configurations separated by 40 
sweeps after a thermalization of 4000 sweeps at (3 V = 0.60, 0.64 and 0.68 on 8 4 lattice. We 



used the DeGrand-Toussaint scheme || for locating monopole currents in the lattice gauge 
field configurations to obtain an ensemble of monopole configurations. The statistical errors 
were estimated with the jackknife method. 

The coupling constants fa ~ fa determined are plotted in Fig.l in comprison with the 
theoretical values, fa agree well with those of the theoretical values, whereas there are small 
discrepancies with respect to f2 ~ fa- The discrepancies come from the truncation of the 
terms of the action taken. Really, if we include more quadratic terms, they disappear. 

Next we applied the method to SU (2) lattice gauge theory. Monte-Carlo simulations were 
done on 12 4 , 14 4 , 16 4 , 18 4 , 20 4 and 24 4 lattices for 2.3 < < 3.0. The U(l) link fields are 
extracted after the abelian projection in the MA gauge. The DeGrand-Toussaint scheme 
@ is used to get monopole configurations. Long-distance behaviors are expected to be 
important in the confinement phase of QCD. Hence we consider extended monopoles of the 
type-2 p0^p2|. The n 3 extended monopole of the type-2 has a total magnetic charge inside 
the n 3 cube and is defined on a sublattice with the spacing b = na, a being the spacing 



of the original lattice [23|. The definition of the type-2 extended monopoles corresponds to 
making a block spin transformation of the monopole currents with the scale factor n. We 
call the sublattice as a renormalized lattice. We derived the coupling constants for l 3 , 2 3 , 3 3 
and 4 3 extended monopoles. 

Our results are summarized as follows: 

l)The coupling constants fa are fixed for not so large value of b ,i.e., b < 1.3 x 10 _2 (A L ) _1 . 
The iteration does not converge in the case with larger b. The coupling constants determined 
are almost independent of the lattice volume as seen from Fig.2. We see fa is dominant and 
the coupling constants decrease rapidly as the distance between the two monopole currents 
increases. 



2)The energy of a monopole loop (with a unit charge |fc M (s)| = 1 ]23]) of length L may 
be approximated by a self-energy part faL, whereas its entropy grows like Lln7 when L is 
large Q . If fa < ln7, the entropy dominates over the energy and the monopole condensation 
occurs. We plot fa versus f3 for various extended monopoles on 24 4 lattice in comparison 



with the entropy value ln7 for the infinite volume in Fig. 3. Each extended monopole has its 
own (3 region where the condition fx < ln7 is satisfied. When the extendedness is bigger, 
larger (3 is included in such a region. 

3)It is interesting to see the relation between the monopole condensation and the de- 
confinement transition in a finite volume. In ref. |25||, Polyakov loops are measured on 



symmetric lattices and the value T c = (25.8 ± 1.1)A^ is given for the critical 'temperature'. 
The critical coupling (3 C is fixied from the value for each lattice volume. For example, (3 C is 
2.60 on 14 4 and 2.81 on 24 4 lattices, respectively. 

Our data indicate the following features irrespective of the original lattice volume. If 
the renormalized lattice is larger than 7 4 , the /i value is above the ln7 line at f3 c . It seems 
to cross the ln7 line just at f3 c when the renormalized lattice is about 7 4 . See Fig. 2 for 2 3 
extended monopoles at f3 — 2.6 which is just (3 C of 14 4 lattice and in which fx takes the 
value about ln7. Also ,in the Villain case shown in Fig.l, the energy-entropy arguments give 
the correct critical coupling on a lattice as small as 8 4 . We may assert the above energy- 
entropy arguments can be used when the renormalized lattice is larger than 7 4 . Then we 
see condensation of some extended monopoles occurs in the confinement phase. 

On the other hand, when the renormalized lattice is smaller than 7 4 , such simple argu- 
ments may not apply as seen from small discrepancy between the ln7 cross point of /i in 
the 4 3 case and (3 C in Fig. 3. There is a possible entropy decreasing effect due to the periodic 
boundary condition on a small lattice. Such entropy calculations are very important to know 
if the confining phase and the monopole-condensed phase are exactly the same as expected 
2§. 

4) The behaviors of the coupling constants are different for different extended monopoles. 
But if we plot them versus h , we get a unique curve as in Fig.4. The coupling constants 
seem to depend only on b, not on the extendedness nor (3. This suggests the existence of the 
continuum limit and the monopole action in the limit may be similar to that given here. 

5) From Fig.4, we see the self-energy per monopole loop length j\ on the renormalized 
lattice decreases as the length scale b increases. A critical length b c ~ 5.2 x lO- 3 ^)- 1 



exists at which the f\ value crosses the ln7 line. The b dependence and the lattice-size 
independence of fi suggest the following picture of the QCD vacuum on the standpoint of 
the monopole condensation. When the physical volume of the system l A (I = Na((3) on iV 4 
lattice) is large enough, that is, I ^> b c , one can take a large renormalized lattice with a 
spacing b > b c . The entropy of monopole loops dominates the energy on such a renormalized 
lattice and the monopole condensation occurs. On the other hand, when the physical volume 
is small, the renormalized lattice must have a small lattice size or a small lattice spacing or 
both. A small lattice size would lead to a possible entropy decreasing and a small spacing 
makes the self-energy large. Thus the monopole condensation does not occur in a small 
physical volume. The transition from the monopole-condensed phase to the normal phase 
takes place on a finite physical volume / 4 . l c must be several times of b c . If it corresponds 
to the deconfining phase transition, we get l c ~ 7b c . The existence of l c means that there 
are always both monopole-condensed and uncondensed phases on a finite (in a lattice unit) 
lattice. 

6) The monopole action may be fitted by 

S = ^mobk^k^s) 

1 Air 

where g{b) is the SU(2) running coupling constant 

9(r2 =2^ ln( ^) + 4i^ lnln (^ < 4) 



The scale parameter determined is A ~ [p7 |. D(s) is a modified lattice Coulomb 



propagator g. For details, see ref. p8fl . The solid line is the prediction given by the action 
with the parameters written in the figure. This form of the action is predicted theoretically 
by Smit and Sijs ||. The existence of the bare monopole mass mo and the running coupling 
constant g(b) is characterisitic of the action in comparison with that of compact QED. 

In summary, QCD confinement may be interpreted as the dual Meissner effect due to 
monopole condensation in the maximally abelian gauge. Our data is the first that shows the 

6 



possible occurence of the monopole condensation in QCD. In the second part of this note 
p9| , it will be shown that the monopoles alone are responsible for the string tension also in 
SU(2) QCD as is in compact QED 0. Details will be published elsewhere |28| . 

We wish to acknowledge Yoshimi Matsubara for useful discussions especially on entropy 
decreasing effects of monopole loops. This work is financially supported by JSPS Grant-in 
Aid for Scientific Research (c) (No. 04640289). 
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TABLES 

TABLE I. The form of monopole action adopted. In addition to the following, all terms in 



which the relation of the two currents 


is equivalent are added. 


i 


interaction 


1 


fci(0,0,0,0) 2 


2 


fci(0,0,0,0)fci(l,0,0,0) 


3 


fci(0,0,0,0)fci(0,l,0,0) 


4 


fci(0,0,0,0)fci(l,l,0,0) 


5 


fci(0,0,0,0)fci(0,l,l,0) 


6 


fci(0,0,0,0)fci(2,0,0,0) 


7 


fci(0,0,0,0)fci(l,l,l,0) 


8 


fci(0,0,0,0)fci(0,l,l,l) 


9 


fci(0,0,0,0)fci(l,l,l,l) 


10 


fci(0,0,0,0)fci(2,l,0,0) 


11 


fci(0,0,0,0)fc 2 (2,l,0,0) 


12 


fci(0,0,0,0)fci(2,l,l,0) 
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FIGURES 

FIG. 1. Coupling constants fi versus (3 in the Villain model. f\ (circle), fi (square), 
(diamond), f± (triangle) and /s (cross) are plotted. Solid lines denote the theoretical values. 

FIG. 2. Coupling constants fi versus lattice size . f\ (circle), f'2 (square), (diamond), 
(triangle) and /s (cross) are plotted. 



FIG. 3. Coupling constants fi versus j3 for 2 3 ,3 3 and 4 3 extended monopoles on 24 4 lattice. 



FIG. 4. Coupling constants f\ versus b 
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